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Abstract 

Some aspects of asymptotic freedom are discussed in the context of a simple two-particle non- 
relativistic confining potential model. In this model asymptotic freedom follows from the similarity 
of the free-particle and bound state radial wave functions at small distances and for the same an- 
gular momentum and the same large energy. This similarity, which can be understood using simple 
quantum mechanical arguments, can be used to show that the exact response function approaches 
that obtained when final state interactions are ignored. A method of calculating corrections to this 
limit is given and explicit examples are given for the case of the harmonic oscillator. 

PACS numbers: 13.60.Hb, 13.60.-r, 25.30.Dh, 25.30.Fj 
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I. INTRODUCTION 



Asymptotic freedom was discovered to be a property of non-Abelian gauge theories nearly 
30 years ago , and has been used to explain scaling in electroproduction p . A version 
of asymptotic freedom also holds trivially in non- relativist ic quantum mechanics since the 
Born series converges more and more rapidly with increasing energy. In particular, final 
state interactions can be ignored in the response functions for inclusive scattering from 
condensed matter systems at large momentum transfer [Q. This result follows easily from 
simple quantum mechanical arguments, although the calculation of corrections can be quite 
complicated |^. The case of confined systems is a bit more subtle, however, since the 
particles in the final state are always in discrete excited states of finite (possibly large) 
size and thus never really free. The case of non-relativistic confinement has been discussed 
extensively, however, as a toy model for asymptotic freedom and quark-hadron duality 
in the structure function for electroproduction, where the production of resonances at low 
energies and momentum transfers averages to the smooth scaling curve observed at large 
momentum transfers. This is now known to occur quite locally, essentially for each resonance. 
An explanation in terms of a relativistic quark model has been given recently 0], while Close 
and Isgur have used a simple non-relativistic model to explain one aspect of duality: how 
the square of a sum in the amplitude for production of resonances becomes approximately the 
sum of squares required for duality. A similar non-relativistic model had been used earlier 
by Greenberg to investigate the scaling limit in terms of the Bjorken scaling variable. 

In this paper we also use a two-body model to discuss non-relativistic asymptotic freedom. 
Such a model provides at best an approximate description of mesons containing heavy quarks 



pHI , p/L| , but it is simple enough that some features of the large momentum transfer limit 



of the response function can be illustrated very concretely. In the context of the model, 
asymptotic freedom follows because of a simple relation between the radial wave functions 
describing bound states in the confining potential and the spherical Bessel functions which 
describe free particles of the same energy and angular momentum. Not only do these wave 
functions have the same shape for high enough energy and small enough separation: their 
relative normalizations are such that a form of duality follows once one properly allows for 
the comparison of a discrete sum with an integral. 

The paper begins with the definition of the response functions of the model, showing how 
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they are related to integrals over the final state radial wave functions. Both the exact case, 
where the final state wave functions are just the bound state wave functions for the discrete 
final states, and the approximate case using free-particle (plane-wave) wave functions are 
considered. The close relation between the two sets of wave functions for the same angular 
momentum, similar large energies, and small radii are then described in Sect. 3 for fairly 
general confining potentials, and the implications of these results for local duality are dis- 
cussed. A method for calculating corrections to this simple relation is also given, leading 
to a discussion of the conditions under which asymptotic freedom is a good approximation. 
The results are applied in Sect. 4 to the case of the harmonic oscillator, where many of the 
expressions become simple analytic functions. The results are summarized in Sect. 5, which 
also contains suggestions for further work in this area . 



II. STRUCTURE FUNCTIONS 



We use the model of 0, consisting of two particles of equal mass m and reduced mass 
/i = m/2 carrying charges Ci and 62- This system is initially in its ground state in the 
spherically symmetric confining potential V{r) and, after being hit by a "scalar photon", 
transferring momentum q to one of the particles, it makes the transition to an excited state 
\nr,l,m > with energy En^^i . The probability for this transition is proportional to the 
response or structure function 

Fn.,i,M = \<nr,l, m\ [eie^'-'i/^ + e2e~'''-^/^]\0, 0,0 >\\ (1) 

If the polarizations m of the degenerate final states of a given and / are not measured 
one needs only the sum 

Fn,.A<l) = E Fr,^,l,M (2) 
m=—l 

Using completeness it is easy to see that 

Ei^n.,/(g)=e? + e2 + 2eie25(g), (3) 

•nr,l 

where 

5(g) =< 0,0, 0|e'"-^|0, 0,0 > (4) 
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is the ground state form factor. Using standard manipulations one can write Fn^^i{q) in terms 
of the square of a radial integral involving the bound state radial wave functions Un^^i{r): 

FnrM = [e? + el + 2eie2(-l)^](2/ + 1) r„,,Kg)', (5) 

where 

rnrM)= drun,,i{r) ii{qr/2)uofl{r), (6) 
Jo 

with ii{x) a spherical Bessel function. The (—1)' in ^ shows that the interference terms will 
tend to cancel when states of adjacent values of / are included in a sum [Q. 

If the interaction of the two particles in the final state is completely ignored, as would 
follow from the assumption of asymptotic freedom, the final states can be labelled by either 
the relative momentum k of the two particles or its magnitude k and the angular momentum 
/, m. Both will be useful below, so we define 

F(k, q) = I < k| [eie''-'^'^ + eae-'^'i/^] |0, 0, > 1^ (7) 

and 

i^^,n^(A:,q) = | < fc, /, m|[eie*^ 'i/2 + e^e'^-'i/^jlO, 0, > \\ (8) 
With these definitions 

F(k, q) = e?02(|k - q/2|) + e^02(|k + q/2|) + 2eie20(|k - q/2|)0(|k + q/2|), (9) 

where is the ground state momentum space wave function which is large only when p 
is less than 1/ro, where tq indicates the size of the ground state wave function, and 

I 

Fi{k, q)= Fi,m{k, q) = [ej + + 2eie2(-l)'] (2/ + 1) n^k, qf, (10) 

m=—l 

with 

ri{k,q)=2 drji{kr)ji{qr/2)uofi{r). (11) 



Here, to emphasize the similarity to the bound state case, we have introduced the Ricatti- 
Bessel functions j{x) = xj{x). Summing Fi{k, q) over / produces the same structure function 
obtained by integrating -F(k, q) over all directions of k : 

oo „ 

F{k, q) = Y. m, Q) = {k/2nf / dQ^F{k, q). (12) 
1=0 ■' 
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The completeness relation now take the form 

/ {dk/27i)F{k,q) = el + el + 2eie2Siq). (13) 



Several properties of the structure functions can be read off from the expressions above. 
In the free case the structure function is large when k ^ ±q/2 , and, for large q, the 
interference term cannot be large. The corresponding result for Fi{k, q) can be seen in Eqn. 



n], since here the two spherical Bessel functions are exactly in phase only if k ^ g/2, so 



that ri{k, q) will be maximum here and begin to decrease significantly when \k — q/2\rQ > 1, 
where tq indicates the size of the ground state wave function. Using Eqn. ^, a simple 
change of variable in the integral over the direction of k shows that, in the large q limit, 
F{k, q) is simply the ground state probability distribution for the component of k in the q 
direction |jl2[ . This means that, in this limit, F{k, q) depends only on the scaling variable 
k — q/2 ^ (k"^ — q^ /^)/qi which is the non-relativistic version of Bjorken scaling . 

In the confinement case the structure function Fn,.,i can be large for large q only if the 
oscillations of the radial wave function Un^^i{r) are in phase with those of the spherical Bessel 
function ji{qr/2), so that asymptotic freedom will hold only if these radial wave functions 
have a form similar to that of j{kr). In the next section this requirement will be studied in 
more detail and it will be shown that the two radial wave functions have almost the same 
shapes for the same energy and angular momentum. 



III. RADIAL WAVE FUNCTIONS 



Assuming V{r) smoothly as r 0, the bound state radial wave function must 
have the same shape as the free particle radial wave functions for the same / and for high 
enough energy and small enough r since they satisfy approximately the same wave equation 
and both vanish at r = 0. In this section we use the WKB approximation to discuss the 
normalization of the bound state radial wave function and show that the result obtained 
guarantees asymptotic freedom. 

In the WKB approximation the radial wave function in the classically allowed region is 
simply 

Un,,i{r) ~ ^nr,i cos[ T dv' kn,,i{r') - 7T / 4] / ^fk^Ar) , (U) 
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where 

KAr) = ^kl,i - 2fiV{r) - {I + 1/2Y/t^ (15) 
is the classical radial wave number at r for a system with energy 

Enr,i = klj2f^, (16) 

with the usual replacement /(/ + 1) — (/ + 1/2)^ [0, and the classical turning points r_ 
and r+ are determined by the condition kn,.,i{r) = 0. The normalization constant A/'„^,/ can 
be estimated by assuming that, for high energy bound states, the average of the square of 
the cosine is close to 1/2 and that the contributions to the normalization integral from the 
classically forbidden regions is negligible. Normalization then requires 

r2 ' 



l^AflJ dr/{2K^,i{r)). (17) 

The integral is clearly proportional to the period of the classical motion, and it is easy 
to show from the Bohr-Sommerfeld condition that it is inversely proportional to the level 
splitting. The final result can be written 

Ml^i ^ 2/x(E„,+i,, - En^,i)/n. (18) 

(Similar arguments can be used to derive the expression for the the local electron density 
used in Thomas- Fermi theory 0].) 

From the discussion above, the shape of Un^^r) must match that of the Ricatti- Bessel 
function ji{kn^^ir) at small r . If we assume this match extends into the classically allowed 
region and to kr ^ /, so that the Ricatti-Bessel function takes on its sinusoidal asymptotic 
form, it must be approximately true that 



i{r) ^Afnr,iJi{knr,ir)/Jkn^j. (19) 



This relation holds only for small enough radius since the local wave number knr,i{i^) ap- 
pearing in the WKB approximation will eventually begin to differ significantly from from 
kn^^i. The radial integrals Eq. (|,|l3) required for calculating the structure functions, how- 
ever, involve only values of r such that the ground state radial wave function Mo,o('") is large. 
As shown in Fig. 1, this can be much smaller than the size of higher energy excited states 
and so it might well be possible that the radial integrals in the bound and free cases are 
essentially identical except for the normalization constant i. 




FIG. 1: Oscillator bound state radial wave function for = 5 and / = 2 (solid) and the free 
particle radial wave function(dashed) for the same energy and angular momentum, together with 
the ground state (n^ = 0, / = 0) radial wave function, plotted against the dimensionless radial 
variable x = ^/ar. The heavy curve shows the effective potential in the oscillator well for this 
angular momentum. 

Using the relation |16| between En^^i and knr,i it is easy to show that 

f kdk/{fiAE), (20) 

rir 

where AE is the energy difference appearing in Eq. |18| for energies related to k by Eq. 
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Then, using (H) and (pi]), 

E^",.,,/(?)^ fdkFi{k,q)/i27r). (21) 

These approximate relations hold for any set of n^, 's provided the integral on the right is 
over the corresponding range of k's, in the sense that the range of energies covered on the 
two sides is the same. This of course is just a statement of asymptotic freedom. Note 
that it holds for groups of states with a single definite angular momentum /, so that the 
correspondence is local in angular momentum as well as energy. 

It has been stated somewhat vaguely above that the above results hold for small enough 
r and high enough energy. It would be useful to have a more quantitative estimate of the 
errors and a scheme for correcting them. One method of calculating corrections has been 
discussed by Rinat et al 0. Another, of more direct application to the approach here, 
using radial wave functions, can be obtained from the integral equation for the radial wave 



function 15 



UnrA^) = ji{knr,ir) - {l/kn„i) / dr'[ji{kn„ir)hi{kn„ir') - ji{kn„ir')ni{kn,,ir)] 

Jo 

2fiV{r')un,.,i{r'), (22) 

where hi{x) is the Ricatti-Bessel function of the 2nd kind and Unr,i{^) is an un-normalized 
version of u which equals jiikn^^r) for very small r. (In this equation u is normalizable only 
if knr,i corresponds to an energy eigenvalue.) This integral equation shows clearly that the 
error in the proportionality between u and ji generally increases with r and the strength of 
V , but decreases as kn^^i increases. It can be used to develop a perturbation expansion for 
the correction, the leading correction being simply the integral term in (^2]) with ^^^^/(r') 
replaced by ji{kn^^ir'). This expansion must, of course, diverge at large r, but it can be 
used to estimate the accuracy of the proportionality ( |T9|) and thus of local duality itself. A 
simple change of variables, for example, shows that for an r" potential the leading correction 
decreases as l//cj^^f . 

IV. HARMONIC OSCILLATOR EXAMPLES 

In the case of the harmonic oscillator potential V{r) = ^Kr'^ there are analytic expressions 
for most of the quantities discussed above. It will be useful to define the usual parameters 
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ijj = \jK/fi and a = fiu, in terms of which = uj{2nr + I + 3/2). The radial wave 

functions are then given by 

Un.Ar) = iV..,«^/V+^Lit^/2(x2)e-^V^ (23) 

where -Z^^^^^ is an associated Laguerre polynomial, the dimensionless variable x = y/a r and 
the normalization constant is 

iVn.,/ = [2nJ/r(n, + / + 3/2)]i/l (24) 

Since in this case kn^^i = ^J2a{2nr + / + 3/2) , the relation between the bound and free 
radial wave functions in terms of dimensionless quantities is 

Wn,,/(r)/a'/^ ^ 2ji{^jAnr + 2l + 3x)/[y/^{Anr + 21 + 3)^/% (25) 

An example of this relation is shown in Fig. 1, which shows that the normalized bound state 
and free particle radial wave functions for the same I and energy are indeed almost identical 
near the the origin, in the region where the ground state radial wave function appearing in 
the radial integrals is large. 

For the oscillator the radial integrals have simple analytic expression: 

rn.,/(g) = v/2"^+VK!(2n, + 2/ + l)!!] (p/2)2-+' e-p'/^ (26) 

and 

ri{k, q) = 2y2/€(7r/a)^/^e-("'+*'')/2 (27) 

where p = q/ {2\/oi) and n = k/ ^fa are dimensionless momenta and ii is a modified spherical 
Bessel function of the first kind. These two functions are plotted in Fig. 2 for the particular 
case of = 5 and / = 2. 

With these expressions for r„,,_i(g) and r;(/c, g) we can use (|]) and ([T0|) to obtain expres- 
sions for ;(g) and Fi{k, q). These can be summed over / to obtain the structure functions 
for transitions to all states at a definite energy: 

= [el + el + 2eM-m {l/^^ivV^Te-^''^ (28) 
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FIG. 2: Radial integrals for = 5 and 1 = 2, computed using bound state (solid) and free 
particle (dashed) radial wave functions for the same energy (k = V2nr- + Z + 3/2 = ^/27 ^ 5.2 ), 
as a function of the dimensionless momentum transfer p = q/{2^/a). Note that both curves peak 
near p = 5.2. 

and 

oo 

F{k,q) = Y.Fi{k,q) 

1=0 

= [el + el] 2y^(/€/p)[e-("-P)' - e^("+P)'] 

+2eie28^f^K^e"P'-''\ (29) 

It should be noted that for large momentum transfers the /-dependence of the terms in 
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the sums above becomes Gaussian: {21 + l)exp{~P /2n) in ( pS]) and {21 + l)exp{—l'^ / up) 
in ). This means that the number of terms which contribute significantly to the sums 
is only a few times \/2n or, equivalently, y/np, respectively, so that the difference between 
the finite sum in ( p8D and the infinite sum in (|29| ) is not significant. This limit on the 
angular momenta of the internal states produced can be easily understood by noting that 
the impact parameters involved cannot be greater than the size of the initial bound state. 
The expression for -F„(g) shows how the contributions from alternate energy levels tend to 
cancel for the 6162 interference term. For F{k,q), on the other hand, the interference term 
is always small for large values of momentum transfer. Furthermore, for large p only the 
exp[—{K, — p)"^] term can be large which, as will become clearer below, gives scaling of the 
structure function. These structure functions, of course, satisfy the sum rules 

J2 Fn{q) = / dk/{2n)F{k, q)=el + el + 2e^e2S{q), (30) 

n=0 ''^ 

where here the form factor S{q) = e~^^ , showing exactly how the interference term vanishes 
as q increases 0. 

In comparing the two structure functions Fn{q) and F{k,q) the fact that in the former 
the energy levels are discrete, while in the latter they form a continuum, must be addressed. 
Here the nth excited discrete level is given a half-width a„ in e, the energy in units of uj. 
This is done by defining 

00 

F{a, e, q) = Fo{q)6{e - eo) + E ^n(g)5a„ (e - e„), (31) 

n=l 

where e„ = n + 3/2, and 5a is a finite-width version of the Dirac delta function S : 

5a{e - e„,) = {a/n)/[{e - e„)' + A- (32) 

In general we expect the widths a„ of the individual resonance peaks to increases with n . 
Furthermore, if 

F(e,g) = (a/(27rfc))F(A;,g), (33) 

where e and k are related by e = k'^/{2a), then the integral over e from zero to infinity 
for both F{a,e,q) and F{e,q) equals the right-hand expression in (|30D. (This equality 
is approximate for F{a, e, q) because the tails of the F{a, e, q) extend into the negative e 
region.) The two F-functions are plotted in Fig. 3 for the case where Ci = 62 = e, for 
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FIG. 3: Structure functions as a function of the dimensionless energy e = n + 3/2 = k"^ /2a for 
dimensionless momentum transfers p = q/{2^/a) equal to 4 and 8. The smooth curves are for the 
free case F{e,q) while the curves with peaks are for the bound case F{a,€,q). In the latter cases 
the sharp energy levels have been given widths which increase from 0.2 to 2 as their energy 
increases. 

which the square bracket in (pS]) is alternately 4e^ and : the contributions from the two 
particles cancel or add coherently for odd and even parity states respectively P] . Note that 
the peaks of the smooth(ed) curves appear approximately at e = and that the widths 
in e increases nearly linearly with p. (The heights of the curves therefore decrease inversely 
with p.) Thus plotting pF versus a scaling variable y = {e — p'^/2)/p /, gives curves of a 
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FIG. 4: Structure functions F times p as a function of the scaling variable y for dimensionless 
momentum transfers p equal to 4 (large peaks), 8 (small peaks), and 32 (smooth curve). The sharp 
energy levels have been given a width which increases from 0.2 to 2 as their energy increases. The 
dashed curve is the Gaussian limit of the scaled free particle pF: (2/y^)e~^^. 

nearly Gaussian shape almost independent of p as p ^ oo, as shown in Figure 4. This scaling 
variable is, aside from a trivial factor, just the y variable of West |T^, and corresponds to 
the component of a parton's momentum in the q direction before the collision. 

We can also calculate the first order correction to the approximations in which the bound 
state wave functions for the oscillator are proportional to spherical Bessel functions. The 
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first two terms in the expansion can be written as 

VurA^) = Vnlli^) + Vnllix), (34) 

where 

v^Zi^)=m (35) 

and 

v^^l^ix) = -[l/{Anr + 21 + 3)2] £ dx'[Ux)ni{x') - ji{x')ni{x)] x'^M^'), (36) 

with X = kn^^ir. An indication of the size of the errors and the accuracy of the leading 
correction above for the case rir = 5 and / = 2 are shown in Fig. 5. It is clear that 
the corrections are small in the region where the ground state wave function is large, but 
including them does improve the accuracy of the expression (0), especially at larger values 
of r. Note, however, that the corrections diverge rapidly at very large r. 



V. CONCLUSION 



In the non-relativistic case asymptotic freedom and scaling are automatic for high en- 
ergy transfers provided the confining potential approaches zero smoothly as r — > 0. This 
result is almost trivial since the system is effectively free at small separations. Only the 
relative normalizations of the terms in the sums for the bound states and the integrals in 
the free case require additional argument from the WKB approximation or, equivalently, 
the correspondence principle. 

As is well known, the resulting structure functions are essentially the same for the confined 
and free cases, except for the resonance-like bumps in the former, and for large momentum 
transfers depends only on the scaling variable y = {t — p'^/2)/p ^ k, — p, as shown in Fig. 3. 

It would be of interest to extend these results to more general potentials and to semi- 
relativistic calculations |Ty]. For example, how would the results be modified for the 
frequently-used Coulomb-plus-linear confining potential []TB|, or for various relativistic ex- 
tensions 1^, |18[? Another possible extension would be to many-body systems, including 



the nucleon and nuclei 12 
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FIG. 5: Un- normalized radial wave functions and corrections for = 5 and / = 2 versus the 
dimensionless radius x. The heavy solid curve is the exact radial wave function, while the lighter 
solid and dashed curve are 20 times the difference between j2 and the exact wave function, and 
the first order estimate of this difference given by (]36|), respectively . The heavy dashed curve is 
the ground state radial wave function. 
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